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ABSTRACT 

Jets from active galactic nuclei (AGN) in the cores of galaxy clusters have the potential to be a major 
contributor to the energy budget of the intracluster medium (ICM) . To study the dependence of the 
interaction between the AGN jets and the ICM on the parameters of the jets themselves, we present 
a parameter survey of two-dimensional (axisymmetric) ideal hydrodynamic models of back-to-back 
jets injected into a cluster atmosphere (with varying Mach numbers and kinetic luminosities). We 
follow the passive evolution of the resulting structures for several times longer than the active lifetime 
of the jet. The simulations fall into roughly two classes, cocoon-bounded and non-cocoon bounded 
sources. We suggest a correspondence between these two classes and the Faranoff-Riley types. We 
find that the cocoon-bounded sources inject significantly more entropy into the core regions of the 
ICM atmosphere, even though the efficiency with which energy is thermalized is independent of the 
morphological class. In all cases, a large fraction (50-80%) of the energy injected by the jet ends up 
as gravitational potential energy due to the expansion of the atmosphere. 

Subject headings: cooling flows - galaxies: active - galaxies: jets — hydrodynamics - X-rays: galaxies: 
clusters 



1. INTRODUCTION 

Jets from powerful Active Galactic Nuclei (AGN) 
play a major role in shaping their environment on 
both the large and small scales. In many clus- 
ters of galaxies, the intracluster medium (ICM) bears 
the marks of repeated episodes of AGN jet activity : 
e.g., bubbles (IFabian et al.ll200Ct IMcNamara et alJl2000t 
Blanton ct al. 200 it I Young et al J 12002ft. ghost bubbles 
(IMcNamara et al.l [20011: iHeinz et all 120021: IChoi et all 



20041 iFabian et all 
20051 ) . shells ([Fabian et"atl l2000ft . 



20001), ripple s (|Fabian et alj 12003. 

and filaments. Be- 
sides affecting the gross morphology of the ICM, AGN 
have the potential to dramatically alter the energy and 
entropy budgets of the ICM. This is especially of in- 
terest in relation to the cooling flow problem of galaxy 
clusters (ICowie fc Binnevlfl977t IFabian fc Nulseni Il977t 
IFabian! 11994ft . In this paper, we attempt to determine 
the relationship between the jet parameters (Mach num- 
ber and kinetic luminosity) and the resulting energetics 
and structures in the ICM. 

X-ray observations show that the ICM radiates 
strongly in the X-ray wavelengths due to thermal 
bremsstrahlung radiation. In the central regions of rich, 
relaxed clusters, the ICM has a cooling time shorter than 
the Hubble time (sometimes as low as a few xlO 8 years). 
With such short cooling times, we would expect to see 
large quantities of cool gas and/or star formation in the 
center of rich clusters. Equivalently, the cD galaxy at the 
center of the cluster would be accreting cooled matter 
from the cooling ICM. The density-squared dependence 
of the emission ensures that, if this simple picture were 
correct, the cooling would increase and eventually reach 
a cooling catastrophe in a finite time. 

Observations strongly suggest that large quantities of 
cooled gas are not being deposited in the central galaxy. 
The observed star formation rates in the central cD 
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galaxies are not sufficient to match the inferred coo l- 
ing rates (|Q'Dea et al.ll2004HHicks fc Mushotzkvll2005ft . 
Furthermore, the mass function of galaxies shows a high- 
mass truncation that argues against the continued ICM 
accretion growth of cD galaxies and requires a feedback 
pro cess significantly mo re efficient than star formation 
(see lBenson et al. (2003)). More directly, high-resolution 
X-ray spectroscopy with the Reflection Grating Spec- 
trometer (RGS) on XMM-Newton shows that the ICM 
cools from the virial temperature to approximately one 
third of the virial temperature but reveals an absence of 
plasma below this temperature. The apparent contra- 
diction between the absence of cool gas and the short 
cooli ng times of t he ICM core is the cooling flow prob- 
lem (|Fabianlll994ft . 

These observational results suggest that some form of 
ICM heating is required to offset the radiative cooling. 
Intermittent activity by a central AGN remains an at- 
tractive solution. A large amount of theoretical work has 
been performed on the possible effects of AGN on cooling 
cluster cores. With ever increasing computer power, 
most of the recent work has focused on h ydrodynamic 
mode l s of the AGN / ICM interaction (IClarke et all 
19971 IChurazov etai] 2001: B riiggen fc Raised 120011 . 



Reynolds et alj l2002t" iBasson fc Alexander! 

Omma & Binney 2004; 



fOmma et alJ 
Dalla Vecchia et "all 
Zanni et all 12005ft . 



2001" 



20041 [Robinson eLaJJ 



have also 



2004; 
begun 



Some groups 

initial explorations of the effects of magne tic fields 
(Ro binson et al.ll2004t iJones fc D e Young 200 51), plasm a 
processes (jRuszkowski et al.ll2004HRevnolds et al.ll2005h . 
feedback prescriptions (|Vernaleo fc Revnoldsl l2006ft. and 



realis tic (cosmological) background motions ( Heinz et al.l 
l200l . 

Collectively, this body of simulation work has allowed 
us to explain many of the observed features of AGN/ICM 
interactions, investigate how AGN induced flows mix 
metals within the ICM core, and study entropy injec- 
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tion and heating of the ICM. However, most of the cur- 
rent sets of simulation are performed for a limited set of 
jet parameters and hence it is unclear how to generalize 
the results to the population of as a whole. In this pa- 
per, we present a moderately large set of (axisymmetric) 
high-resolution hydrodynamic simulations of jet/ICM in- 
teractions that survey a wide range of jet powers and jet 
velocities. We study the morphology of the radio-galaxy 
as well as the injection of energy and entropy into the 
ICM as a function of the jet properties. In Section [51 we 
will discuss the details of our setup and the code used. In 
Section 02 we will present the analysis of our simulations 
and their classifications. In Section H] we will discuss the 
results, followed by conclusions in Section [5] 

2. SIMULATION SETUP 

Our goal is to model a relaxed cluster and its interac- 
tion with a central radio galaxy that produces back-to- 
back jets. To do this, we start out with an initial cluster 
that is both spherically symmetric and isothermal, with 
the (adiabatic) sound speed c s = 1 everywhere. The ini- 
tial cluster gas is assumed to follow a /3-model density 
profile, 

= [l + (J-)2]3/4> W 

where the core radius is r$ — 2.0 in code units. 
A static gravitational potential, 

* = -ln(p), (2) 

7 

is set to make the gas initially in hydrostatic equilibrium 
and is assumed to remain fixed in time. This is equivalent 
to a potential set entirely by stationary dark matter that 
dominates the system, so the self gravity of the gas is 
ignored. The intracluster gas, while containing a large 
fraction of the baryons in the cluster, is indeed not a 
significant contributor to the overall cluster mass. 

We use spherical polar coordinates (r, 9, <fi) and impose 
symmetry in the <p direction (d/d<p = 0). The radial 
coordinate varied between 0.1 and 30.0 in scaled code 
units. The comparison between code units and physi- 
cal units is given in Section 12.21 The values of 9 were 
allowed to vary between and 7r (i.e., we are modeling 
the full sphere). Since we have the full range of angles, 
this allows us to have back-to-back jets. Thus we can 
model interaction between the backflows from the jets 
rather than the commonly employed technique of impos- 
ing reflect ing boundary co n dition s on the 9 = ir/2 plane 
(e.g., see ISternberg et al.l (|2007f k IVernaleo fc Reynolds! 
(2006)). We included second order differences in grid 
spacing in the 9 direction for each hemisphere, giving us 
two realizations of the jet/cocoon structure per simula- 
tions. In no cases were there major differences between 
the two sides, confirming that the details of the griding 
do not affect our results. A small circle around r = 
was excluded from the computational domain to avoid 
the coordinate singularity. 

All simulations were run with a n r x rig = 1200 x 600 
grid. In the radial direction, a ratioed grid was used, 
with successive cells increasing by a factor of 1.003. The 
9 direction also used a ratioed grid increasing by a factor 
of 0.998 from to 0.5-7T and then decreasing by a factor of 
1.002 from 0.57T to n. This provides the greatest number 



of zones near the center and near the two jet axes while 
providing the second-order differences in the 9 direction 
mentioned above. 

Clearly jets are three dimensional structures. How- 
ever, there is a long history of two-dimensional mod- 
eling of jets w hich supports the usefu l ness of two di- 
mens i onal jets (|Carvalho fc O'Dca 2002; Churazov et al.l 
120011: iMizuta et all 120011 : iRevnolds et al.ll2002t) . There 
are however several aspects of the dynamics that can- 
not be captured in two dimensions. There is no way to 
have realistic random motions in the background. This 
lack of random motions produces lobes which are far 
more regular and symmetric than those found in any real 
source. The two dimensional assumption also enhances 
mixing between the high entropy material and ambi- 
ent material. This is because non- axisymmetric Kelvin- 
Helmholtz (KH) and Rayleigh- Taylor (RT) modes are 
clearly not possible in an axisymmetric system. Without 
non-axisymmetric modes, all of the mixing must be done 
by the axisymmetric modes, which are more efficient at 
mixing (|Revnolds et ah 2002). Fortunately, this will be 
most problematic only at late times when there has been 
enough time for significant mixing to occur. 

The hydrodynamic evolution of the jet/ICM system 
was followed using the ZEUS-3D code ([Stone fc Normanl 
[1992 a bl. ZEUS is a fixed grid Eulerian code, explicit 
in time and second order accurate (when using van Leer 
advection as we do). Artificial viscosity is used to handle 
shocks. All simulations were run as ideal hydrodynamic 
cases, neglecting magnetic fields and self gravity. ZEUS 
solves the standard equations of hydrodynamics, 



Dp 
Dt 



pV ■ v = 0, 



Dv 

p— = -VP - pV$, 



where 
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Dt 
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dt 



(3) 



(4) 



(5) 



(G) 



Radiative cooling, represented by the A-term in Equa- 
tion O is neglected in order to concentrate on the hydro- 
dynamic evolution and the thermodynamic response of 
the ICM to the jets. 

2.1. Modeling the jets 

Two, back-to-back jets were injected into the center 
of the model cluster. This was achieved with an inflow 
boundary condition on the inner radial boundary. All 
the jets were highly supersonic with respect to the am- 
bient sound speed and moderately to highly supersonic 
with respect to their internal sound speeds. Direct injec- 
tion allows us to keep the details of jet acceleration and 
central engines outside of the grid, so we are able to solve 
the equations of hydrodynamics consistently within our 
computational grid. 

The injected jets arc initially conical with half-opening 
angles of 15°, and are injected in pressure balance with 
the ICM core. The drop of internal pressure associated 
with subsequent jet expansion leads to external, pressure 
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TABLE 1 

Parameters defining the 26 simulations presented in this 

paper. Also given are the (visual) morphological 
classification of the resulting structure, as judged at 
time t = 2.0. See text for more details. 



Run 


Jet Active Time 




JV4j e t 


Morphology 


1 


1.0 


26.47 


10.55 


Cocoon 


2 


1.0 


52.94 


21.1 


Cocoon 


3 


1.0 


26.43 


15.89 


Cocoon 


1 


1.0 


26.31 


7.39 


Cocoon 


5 


1.0 


68.97 


10.46 


Cocoon 


6 


1.0 


23.32 


10.54 


Cocoon 


7 


1.0 


52.94 


10.55 


Cocoon 




0.5 


26.47 


10.55 


Cocoon 


g 


1.5 


26.47 


10.55 


V_' ULUUll 


10 


1.0 


8.91 


6.47 


Cocoon 


11 


1.0 


17.91 


9.48 


ULUUll 


12 


1.0 


19.98 


15.0 


Cocoon 


1 3 
xo 


1 n 

1 A) 




20 02 


Cocoon 


14 


1.0 


13.24 


21.1 


Non-Cocoon 


15 


1.0 


13.23 


10.55 


Non-Cocoon 


16 


1.0 


26.47 


21.1 


Non-Cocoon 


17 


1.0 


6.62 


10.55 


Non-Cocoon 


18 


1.0 


9.96 


8.99 


Non-Cocoon 


19 


1.0 


21.01 


13.0 


Non-Cocoon 


20 


1.0 


30.1 


18.02 


Non-Cocoon 


21 


1.0 


24.98 


19.0 


Non-Cocoon 


22 


1.0 


52.99 


5.28 


Unresolved 


23 


1.0 


26.44 


3.71 


Unresolved 


24 


1.0 


26.43 


5.27 


Unresolved 


25 


1.0 


10.05 


5.01 


Unresolved 


26 


1.0 


12.15 


6.54 


Unresolved 



driven recollimation and internal shocks, ultimately re- 
sulting in some new, self-consistently determined opening 
angle. Given these (fixed) assumptions, the properties of 
a given jet are parameterized by the kinetic luminosity 
Lkin and internal Mach number of the injected jet mat- 
ter. These are related to the injection density p, pressure 
p and velocity v by 

M jet = V J- = Vj ^=, (7) 

Lkin = \pj v ^ r in- ( 8 ) 

where r- m is the radius of the "injection nozzle" for the 
jet. The jet parameters for the base run (Xkm = 26.47 
and M- iCt — 10.55; pj = 0.0 1 and Vj = 105 . 5) ar e the 
same as the run presented in lRevnolds et all (|2002f ). 

At total of 26 models were run, where only the kinetic 
luminosity and/or the Mach number of the jets were var- 
ied. In two of these simulations, the jet lifetimes were 
also varied. The full list of simulations is given in Ta- 
ble!]] Internal Mach numbers in the range 3.7-21.1, and 
kinetic luminosity in the range 6.6 — 68.8 were explored. 
The choice of parameters used was based on raising or 
lowering the Mach number and/or kinetic luminosity by 
a factor of two from one of the previous runs. 

2.2. Comparison to Real Units 

Throughout this paper, quantities are quoted in di- 
mensionless code units unless explicitly stated other- 
wise. Since we are solving the equations of ideal hy- 
drodynamics (Equations [3][ni) with no additional physics 
added, there is no unique set of physical scales to our 



problem. Each simulation may therefore be compared 
to several different sets of physical scales. Following 
Reynolds et al. (2002), we shall quote two scalings for our 
simulations (although there is actually an entire three- 
parameter family of scalings). For our rich cluster scal- 
ing, we set the core radius ro = 100 kpc, meaning one 
code unit in r is equal to 50 kpc. Such clusters are hot, 
with sound speed, c s = 1000 km s -1 . This, along with 
the typical central number density no = 0.01 cm -1 , gives 
a time unit of 50Myrs and a kinetic luminosity unit 
of Lkm = 3.5 x 10 44 ergs _1 . We can also consider a 
poor cluster scaling, r$ = 50 kpc, c s — 500 km s" 1 , and 
no = 0.1 cm -1 . These give a length unit of 25 kpc, a 
time unit of 10 Myrs and a kinetic luminosity unit of 
L km = 4.4 x 10 42 ergs" 1 . 

3. ANALYSIS AND RESULTS 

The primary diagnostic for distinguishing injected jet 
material from ambient and shocked-ambient material is 
the specific entropy index, 

P 

o=—. (9) 

pi v ; 

The ambient material (both undisturbed and shock 
heated) has a low specific entropy index (which ranges 
from 0.6 to 9) and the injected jet material has a specific 
entropy index on the order of 1000 for our base simulation 
(run 1). The specific entropy of injected material does 
vary with jet power, going from near 1 for the weakest 
jets to 10 7 for the most powerful. With only two ways to 
change entropy (increase with shocks and increase or de- 
crease through numerical mixing) , specific entropy gener- 
ally provides a powerful way to distinguish injected (and 
shocked) material from the background gas. A cutoff of 
a = 10 was set to separate jet material from ambient ma- 
terial. This is close to the highest ambient value so we 
are not likely to artificially miss any material that should 
count as part of the cocoon/high entropy material. We 
use the same entropy-index cutoff for all simulations to 
allow for direct comparisons to be made. This value was 
chosen to avoid counting ambient material which has the 
same range of values in all simulations. In the case of 
the weakest jets, this entropy index cutoff cannot pick 
out the pre-shocked jet material. Once the jet material 
passes through the shock terminal shock, however, it is 
always clearly distinguished from ambient material by a 
a = 10 cutoff. For a further discus s ion of the entropy 
cutoff diagnostic, see lRevnolds et ail (|2002ft . 

3.1. Morphological Classification 

Based on the entropy maps, the simulations were vi- 
sually placed into three classes. Simulations with a dis- 
tinct, well-defined monolithic cocoon of jet-originating 
plasma were labeled as "Cocoon bounded sources" (or 
simply "cocoon sources"). We consider a cocoon to be 
a jet inflated region of high entropy material with a co- 
herent structure that is well separated from the low en- 
tropy background and relatively safe from mixing (un- 
til late times). Simulations which showed faint, wispy 
plumes of jet-originating material instead of well defined 
cocoons were labeled "Non-Cocoon bounded sources" 
("non-cocoon sources"). The final group were those sim- 
ulations that do not appear sufficiently resolved for an 
accurate classification. In the discussion below, we shall 
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abbreviate the term "jet-originating material" to "jet 
material" . 

An example of a typical cocoon source can be seen in 
Figure Q] Figure [2] shows a typical non-cocoon source. A 
poorly resolved simulation can be seen in Figure [3] We 
also observe that it is possible to separate the cocoon- 
bounded from the non-cocoon bounded sources by use of 
density contours. Figures [4] and \5\ show density contours 
for a cocoon bounded and non-cocoon bounded source 
respectively. Simulations in the cocoon category appear 
to have the high entropy jet material in a distinct region 
separated from the background by a solid density contour 
and remain that way until fairly late times. Non-cocoon 
sources do not have such distinctions and the jet material 
seems mixed in with the background even early on. 

These classifications were performed at the same time 
for all simulations, t = 2.0. At this time, the structures 
have evolved passively for the same amount of time that 
they were driven by the jet. However, our classification 
does not appear to depend sensitively on the choice of 
time, although there are some limitations. It the very 
earliest stages of the active jet phase (t < 0.2 — 0.4), the 
jets are still pushing out from the inner boundary of the 
simulation (i.e., the morphology is strongly determined 
by the artificial aspects of the simulation) and classifi- 
cation is not meaningful. At late times (t > 10 — 15), 
mixing has disrupted all of our simulated sources and, 
again, eliminated any morphological distinction. 

Using this classification, we can study how morphol- 
ogy depends on L- ]et and A4j C t- Our results are shown in 
Fig. O We see three separate regions on this parameter 
space. At the high kinetic luminosities and lower Mach 
numbers, the runs are unresolved. From this point on, 
the unresolved sources were excluded from all analysis. 
Low luminosity sources seem to fall in the non-cocoon 
category, almost regardless of Mach number (except that 
at the high Mach number they can have a higher ki- 
netic luminosity) . The middle region contains the cocoon 
bounded sources. In Fig. [3 we show the parameter space 
again, on a density-velocity plot. Although the overlap 
between the morphologies is more apparent in this pa- 
rameter space, we can see that there appears to be a 
continuum of morphologies with the fastest and light- 
est jet producing unresolved sources at one end and the 
slow heavy jets producing non-cocoon sources. Again, 
the cocoon sources occupy the middle ground between 
the other two. A discussion of the physical origin of the 
cocoon will be reserved for the Discussion (Section |4|). 

As revealed by Fig. [6l the line separating co- 
coon bounded and non-cocoon bounded sources on the 



(Lj, M.j e t) is approximately linear (i.e., Lj 



DC 



M jot ). 



This can be understo od from the analyti c mod els of co- 
coon expa nsion from | Begelman k Cioffil (fl989h . Equa- 
tion (4) of lBegelman fc Cioffil (llQSflj T states the the pres- 
sure of a cocoon is given by 



Pc 



(L^VjPaAh) 1 / 2 



(10) 



where Ah is the area of the working surface over which 
the jet deposits its momentum, A c is the cross-sectional 
area of the cocoon, and p a is the density of the ambient 
ISM/ICM. Well defined cocoons can only be sustained if 
p c exceeds the pressure of the ambient medium, p a . Not- 
ing that the axial ratio of the cocoons are approximately 



constant (i.e., Ah/A c ~ constant), and that the density 
and pressure of our model ambient ICM atmosphere out- 
side of the ICM core in the same in all simulations and 
drops off as r -3 / 2 , the condition that p c > p a can be 
written as the condition 



(ikin^) 1/2 > T c 



(11) 



where T cr j t is a weak function of r (T oc r 1 / 4 ) and hence 
time. Since the jets are injected with fixed initial pres- 
sure, it can be shown that initial jet velocity depends 
on the kinetic luminosity and initial internal Mach num- 
ber according to Vj oc Lkin/A / l? t . Thus, we can use 
can. Illl to see that the line separating cocoon from non- 
cocoon sources (i.e., where the cocoon just comes into 
pressure balance with the ambient ICM) can be written 
as Ljet oc TA^jot- 

3.2. Mass Distribution and Energetics 

To examine the mass distribution and energetics dur- 
ing the simulations, several quantities were calculated at 
each output time. The total mass of jet material at each 
time is given by: 



M, 



cocoon 



P dV 



(12) 



where C is the region with a > 10. The internal, kinetic, 
and potential energies were calculated for both the jet 
and the ambient material, 



E 
E i 
El 



C,A 



C,A 



C,A 
pot 



= ^lJc,A PdV > 
= I fc,A P v2dV > 
= -J C ,AP* dV - 



(13) 
(14) 
(15) 



where A is the region with a < 10. In all cases, the 
integral is taken over only the region of interest (i.e., 
cocoon or ambient only, as determined by the entropy 
cutoff). We then subtract the initial energies in order to 

derive the change of energy AEf^ t , A£^' n , and A_E p0 ' 4 . 

The results of these calculations for a representative co- 
coon bounded simulation are shown in Figure[5] FigureJH] 
shows similar plots for a non-cocoon bou nded simulation. 
These figures are similar to Figure 4 of iRevnolds et all 
( 2002) and primarily differ due to our increased outer ra- 
dius (and the plots are continued until later times). The 
upper-left panel in each figure shows the total mass of 
the cocoon material compared to the amount of injected 
material. In both cases, this increases until the jet turns 
off and then begins to decrease. For the cocoon-bounded 
source, this decrease is mostly steady for the entire life- 
time of the simulation. For the non-cocoon source, the 
drop is much greater early on, so that about halfway 
through the simulation most of the high entropy mate- 
rial has mixed with the background and thus the de- 
cay rate has lessened. This enhanced mixing is a direct 
consequence of the lack of a well defined contact dis- 
continuity in the non-cocoon sources. Note that the co- 
co on mass does not corr espond to the cocoon mass plot 
in Reynolds ct al. ( 2002) as we arc interested in the total 
cocoon mass, not the change in cocoon mass. 

Also shown in Figs. [5] and [5] (upper-right panel) is the 
length of the region containing jet material. In the co- 
coon bounded sources, the length of the cocoon evolves 
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Entropy map of central regions of typical cocoon bounded source (run 1) at t=2.0. 
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Fig. 2. — Entropy map of central regions of typical non-cocoon bounded source (run 21) t=2.0. 




Fig. 3.— 



Entropy map of central regions of typical unresolved source (run 23) at t=2.0. 
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Fig. 4. — Density contours for inner half of cocoon bounded source (run 1) at t=2.0 




Fig. 5. — Density contours for inner half of non-cocoon bounded source (run 21) at t=2.0 
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in a smooth manner, gradually decelerating as the co- 
coon comes into pressure balance with the ambient ICM. 
The discontinuity at late times occurs well into the phase 
in which the cocoon has transformed into two buoy- 
antly rising plumes and corresponds to the complete mix- 
ing/dispersal of the leading part of the plume. The evo- 
lution of the length of the non-cocoon bounded sources, 
on the other hand, suffered a sharp break at the time 
that the jet is shut off. 

The two lower panels of Figs. [H] and [5] show the change 
in energies for the jet material and the ICM material re- 
spectively. For the jet material, there is little difference 
in the time-dependence of the energies for the two classes 
with one exception. In the non-cocoon case, the changes 
in internal and kinetic energy peak at a comparable value 
although the kinetic drops to almost zero very soon af- 
ter. For the cocoon case, the change in kinetic energy 
never reaches a comparable level to the change in inter- 
nal energy. Physically, this is close to the heart of the 
difference between the cocoon and non-cocoon sources. 
In each set amount of energy is injected by the 

jet (along with a set amount of mass as seen in the poten- 
tial change which is not very dramatic). The evolution 
of the system then determines how this energy is split 
between the internal energy and the kinetic energy (and 
eventually how much is transfered to the ICM). In the 
case of the cocoon, the hot jet and shocked material is 
kept separate from the background, and goes almost en- 
tirely into the internal energy of the shocked jet gas. For 
the non-cocoon case, the energy is split nearly evenly 
between the internal and the kinetic energy. This means 
there is no longer enough internal energy available to in- 
flate a cocoon of shocked gas; more of the energy goes 
to the kinetic energy of the mostly bare jet (with wispy 
areas of hot, shocked gas ar ound it). This is c onsist ent 
with the analytic estimate of ICioffi fc B londin (1992). 

The time-dependence of the ICM energies show a 
greater difference between the two cases. The potential 
energy suffers a greater change (within a few times the jet 
lifetime) in the cocoon case due to the well-defined shell 
of ICM that is lifted up by the expanding cocoon, while 
the non-cocoon case has a much more gradual change 
in the potential energy. The internal energy in the co- 
coon case drops back down to the initial value around the 



same time the potential peaks while the non-cocoon case 
has what looks like a slow trade-off between the two. Fi- 
nally, the kinetic energy changes are comparable in the 
two cases, as there appears to be a similar amount of 
disturbance in the ICM regardless of the nature of the 
inflated structure. 

In its standard configuration, ZEUS is not a strictly 
energy conserving code. We do not however believe that 
this is a major impact on our results. Fig. 1101 shows that 
during the majority of time during the simulation, total 
energy is conserved to a sufficient degree. At very early 
times, the jet injects hot material onto the grid which 
is clearly not conservative. At very late times, mate- 
rial (and sound waves carrying energy) leave the outer 
boundary of the grid. During the time in between, the 
energy change remains flat with only minor variations (at 
the 10 -5 level), showing that energy is mostly conserved. 

3.3. Entropy evolution and thermalization efficiencies 

We now examine the effect of the jets on the entropy of 
the background gas. First, we define the specific entropy, 

AS = S X -S =logai -logao, (16) 

where a is the specific entropy index from Equation [9l 

Figure [Tl] shows the entropy difference (compared to 
the initial condition) at the final output of the simula- 
tions. At this point, the systems were allowed to evolve 
passively for nineteen times the active lifetime of the jet. 
For each simulation, the total entropy at each radius was 
calculated and then averaged over the angular coordi- 
nate. The split between the different morphologies is ev- 
ident in this figure. Within the core, the cocoon sources 
have a higher AS overall. Further discussion of this issue 
will be reserved for Section [U 

As a final look at the effect of the jets on the clus- 
ter energetics, we show the efficiencies for conversion of 
the injected energy into various forms within the ICM as 
measured at the final time for both the cocoon and non- 
cocoon sources. This is defined as the change in each 
type of energy for the ICM material divided by the total 
amount of energy injected by the jet at the final out- 
put time of the simulation. Figures [121 El an d El show 
the efficiency for conversion of injected energy into ICM 
internal energy, kinetic energy and potential energy, re- 
spectively. There is no clear segregation between the two 
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Fig. 8. — Mass and Energy vs. time for cocoon bounded simulation (run 1). Starting from the top left and moving clockwise, the panels 
are: Mass of the cocoon material divided by the injected mass, the maximum radial extent of the high entropy material, the change in 
internal, kinetic, and potential energies for the ambient, low entropy material compared to the initial state, and the same energy changes 
for the cocoon material. 
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Fig. 9. — Mass and Energy vs. time for non-cocoon bounded simulation (run 21). Starting from the top left and moving clockwise, the 
panels are: Mass of the cocoon material divided by the injected mass, the maximum radial extent of the high entropy material, the change 
in internal, kinetic, and potential energies for the ambient, low entropy material compared to the initial state, and the same energy changes 
for the cocoon material. 



types although the cocoon sources do appear to cluster 
together somewhat. In all large fraction of the 

injected energy (50-80%) ends as gravitational potential 
energy of the ICM as it expands in response to the AGN 
heating. 

4. DISCUSSION 

In the (resolved) simulations, we see a split in the mor- 
phology of the jet-created structures. As we saw in the 
parameter space plot (Figure O, the presence or absence 
of a cocoon is not related to the jet speed or power in 
a simple matter. So the forms that we see do not tell 
us directly about the parameters of the jet that formed 
them. What they do tell us is something about the way 
the jet interacted with the ICM. 



The cocoon itself is formed by shocks. The supersonic 
jet initially flows freely through the ICM. After traveling 
some distance (usually short on the scale of these sim- 
ulations), the jet will build up enough material in front 
of it that it can no longer flow freely. The material at 
the head of the jet is shock heated and expands. Up to 
this point, the evolution is the same for both cocoon and 
non-cocoon sources. As this shocked material expands, 
the evolution diverges. 

The hot, shocked material expands into the surround- 
ing medium. At the boundaries of the shocked mate- 
rial, fluid instabilities (i.e., the Rayleigh- Taylor (RT) and 
Kelvin-Hclmholtz (KH) instabilities) work to shred the 
contact discontinuity and destroy the forming bubble. In 
the cocoon bounded sources, the bubble is inflated be- 
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Fig. 10. — Fractional change in total energy within the simulated domain as a function of time. All forms of energy (kinetic, internal 
and gravitational) for all of the gas in the computational domain are considered. Shown here is the change in total energy compared with 
the initial time AE divided by the initial total energy Eq. 
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fore it is shredded and continues to inflate until it reaches 
pressure equilibrium with the background. While active, 
the jet deposits energy only at point where it comes into 
contact with the shocked ICM and freshly shocked ma- 
terial flows back from there, continuing to inflate the 
bubble. Once the jet ceases, the bubble floats buoyantly 
from the cluster center. 

In the case of non-cocoon sources, the jet is unable 
to inflate a bubble. If the jet cannot deposit enough 
energy at the contact point to inflate the bubble, the 
fluid instabilities will mix the shocked material with the 
background material, preventing the bubble from ever 
growing. This leaves areas of hot, shocked gas that can 
rise, but never produces a single coherent structure. 

The clearest distinction in the thermodynamics of co- 
coon and non-cocoon sources can be seen in the spatial 
distribution of the entropy enhancement. We find that 
in the core regions, cocoon sources cause a greater en- 
hancement of the ICM entropy than non-cocoon sources 
although there is some overlap. Outside the core (around 
r = 5), both appear the same and cause a slight negative 
change in the entropy. This negative change corresponds 
to lower entropy material that has been pushed upwards 
by the expansion of the core. For the non-cocoon sources, 
this change persists to the edge of the simulation. Some, 
but not all, of the cocoon sources have a spiky, posi- 
tive change when approaching the edge of the grid - 
this corresponds to the remnants of the buoyant cocoon 
that, in some cases, remains intact until the end of the 
simulation. It is important to note that the comparison 
of overall entropy change was done at the final output 
(t = 20). This is long after our classifications were done. 
At this point, mixing has disturbed the structure enough 
that it would no longer be possible to separate the sim- 
ulations into cocoon or non-cocoon categories. The fact 
that we still see a difference at this late a time shows 
that the jets have long lasting effects on the energetics 
of the ICM and that these effects vary based on the jet 
properties in a way that is more complicated than simply 
depending on jet power. 

The energy efficiencies (Figures Q21 [T3l and fT4]) are 
initially puzzling. The non-cocoon simulations generally 
seem to coincide with the highest efficiency cocoon sim- 
ulations for internal and potential energy efficiency. For 
the kinetic energy (Figure [75]). the non-cocoon runs seem 
to have both close to the highest values and the lowest 
values. This is in contrast with the situation for entropy 
where the cocoon simulations have a greater impact. It is 
fairly clear why the cocoon sources have a great entropy 
enhancement overall: shocks. The cocoon is formed by 
shocked gas, so jets that are capable of forming cocoons 
are also generally the ones that shock more ICM (al- 
though this is not the entire story as the rate of shocked 
gas vs. mixed gas matters as well). It is not as im- 
mediately obvious why jets that are good at increasing 
the entropy are not always as good at increasing the en- 
ergy and why there can be a spread in efficiencies for 
jets with similar parameters. To see why, we must go 
back to Equation [51 the definition of the specific entropy 
index (keeping in mind that the pressure, P, is directly 
proportional to the internal energy. For a given change 
in internal energy, the change in entropy can be higher 
(or lower) by having it take place at a lower (or higher) 
density. So if the shocks happen higher in the potential 



well (and hence at a lower density), it is possible to have 
a greater change in entropy for the same amount of en- 
ergy input. The spread in kinetic energies shows that 
it is possible to stir up the gas significantly with both 
the cocoon and the non-cocoon producing jets. Also, as 
stated in the discussion, the cocoon producing jets leave 
more energy available as internal energy which can go to- 
wards shocks while the non-cocoon producing jets have 
more of their energy in kinetic energy which does not go 
to shocks and therefore is not available to increase the 
entropy. It is also the case that different jets will cause a 
different amount of material and energy to be lost from 
the system at the outer edge of the grid. We consider 
it to be reasonable to count this as lost energy as any- 
thing that has not managed to heat the gas by the time 
it reaches the outer edge is clearly not going to be able 
to change the temperature of the system. 

This disconnect between entropy increase and 
energy increase also fit s in with the results 
from IVernaleo fc R eynolds (2006), which showed 
that even with large amounts of energy available, a 
jet can produce large, well formed cocoons yet fail to 
significantly heat the ICM. The large cocoons indicate 
a large change in entropy, but the the central regions 
were still not heated significantly. This helps to reaffirm 
the point that in pure hydrodynamics, it is hard to 
energetically couple powerful jets to the core regions of 
the ICM (where it is needed to solve the cooling flow 
problem!). 

Along with the energetic issues, we may also com- 
pare these simulations to observations in a broad sense. 
Classically, radio galaxies are split into two categories: 
Faranoff-Ri ley Type I (FR I) and Faranoff-Riley Type 
II (FR II) (jFanaroff fc Rilevli 19741 ). Lacking any infor- 
mation on luminosities and radio emission, we cannot 
directly compare our simulations to these classes of real 
radio galaxies, but we can note some interesting corre- 
lations. FR II galaxies encompass the large, back-to- 
back classical double sided sources. Our cocoon-bounded 
sources are reminiscent of the classical doubles. The less 
lobe dominated and often bent FR I sources share a simi- 
lar structure with our non-cocoon sources. Although our 
classifications were done based on the passive phase of 
the system, we can also compare the jets from the active 
phase to the FR I/FR II divide. As seen in Fig. [5] and [71 
slow, heavy, and lower luminosity jets (the non-cocoon 
cases) should correspond to the FR I sources while fast, 
light, and more luminous jets should correspond to the 
FR II type. As FR II gala xies have higher luminosity 
jets (|Ledlow fc Owenl [19961 . this is at least consistent 
with our simulations. Even without a definite correspon- 
dence between the two, we do show that it is possible 
to produce qualitatively different structures with differ- 
ences in the jet parameters without any environmental 
changes. 

Two additional simulations (run 8 and run 9) were per- 
formed where all the parameters were kept the same as 
our standard case (run 1), but the active time of the jet 
was varied. In run 8, the jet had an active time of only 
tj = 0.5 while run 9 had an active time of tj = 1.5. 
In both cases, these simulations fall in the same cocoon- 
bounded category as run 1. This give us some confidence 
that our results are not simply based on the amount of 
total energy inject by the jet, but instead depend on the 
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detailed hydro-dynamic interactions between the jets and 
the ambient material. 

5. CONCLUSION 

We have performed a large number of axisymmetric 
simulations of AGN jets in cluster atmosphere. By lim- 
iting ourselves to pure hydrodynamics with no radiative 
cooling, we are able to carefully study the interaction of 
the atmosphere with the jets and the formation of co- 
coons and bubbles in the cluster. With the continued 
interest in the probable role of jets and bubbles in the 
cooling flow problem, it is important that we understand 
their dynamics as clearly as possible. 

Our main results are summarized in Figures [5] and 111! 
The morphology vs. jet parameters plot shows that we 
are able to produce two distinct classes of structures, co- 
coon bound and non-cocoon bound, by varying the initial 
jet parameters. This distinction does not change with a 
change in total injected energy (as seen when we vary 
the jet lifetime) and is a function of both jet parameters 



(Mach number and kinetic luminosity). We can also vi- 
sually draw comparisons between our two classes and the 
FR I and FR II distinction amongst real radio galaxies. 
The angle averaged entropy change vs. radius shows that 
along with the split in morphology, the effect on the en- 
ergetics of the ICM also depends on the jet parameters. 
We also see that even short lived single burst jets are 
capable of a long-lived enhancement of the entropy of a 
cluster core. 
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